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1.  INTRODUCTION  AND  SUMMARY 


This  is  the  final  rqx)rt  on  ONR  contract  no.  N00014-89-K-0125  for  the  Center  for  Research 
in  Electro-Optics  and  Lasers  (CREOL)  at  the  University  of  Central  Florida  (UCF)  to  provide  basic 
scientific  support  to  the  Innovative  Science  and  Technology  Experimental  Facility  (ISTEF)  at 
Kennedy  Space  Center  (KSC).  The  purpose  of  this  program  was  to  provide  support  to  ISTEF  in 
developing  and  testing  innovative  laser  radar  techniques  and  instrumentation  for  the  discrimination 
and  tracking  of  boosters,  PBV’s,  and  mid-course  targets.  This  reporting  period  covers  the  full 
period  of  the  contract,  from  .^ril  1,  1989  to  December  31,  1990.  Under  this  contract,  specific 
CREOL  support  to  ISTEF  included  the  following  four  tasks: 

Task  I.  One  Micron  Coherent  Laser  Radar, 

Task  n.  Atmospheric  Seeing  Measurements, 

Task  in.  Atmospheric  Induced  Enhanced  Backscatter  Characterization, 

Task  rv.  IR  Focal  Plane  Array  Testing. 

1.1  Task  I.  One  Micron  Coherent  Laser  Radar 

The  initial  goal  of  the  one  micron  task  was  to  mechanically  harden  an  existing  coherent  laser 
radar  system  that  was  to  be  supplied  as  government  furnished  equipment  (GFE).  Accomplishments 
toward  this  goal  were  inhibited  by  the  lack  of  the  delivery  of  the  government  laser  which  resulted 
in  the  adoption  of  an  alternate  program  to  design  and  construct  an  alternate  laser  radar  system. 
Further  delays  were  incurred  as  a  result  of  the  late  delivery  of  a  GFE  transmitter  which  was  to  be 
utilized  in  this  alternate  ladar.  As  a  result  of  these  difOculties  accomplishments  toward  the 
development  of  a  one  micron  coherent  laser  radar  were  hindered.  CREOL’s  accomplishments 
included  the  procurement,  fabrication,  test  and  calibration  of  the  ladar  subsystems  to  include:  1)  a 
local  oscillator,  2)  an  off  axis  Newtonian  telescope,  3)  a  2  Ghz  InGaAs  p-i-n  detector,  4)  a  low  noise 
high  gain  high  bandwidth  preamplifier,  and  5)  assorted  high  power  AR  coated  1.06  optics. 
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12  Task  n.  Atmospheric  Seeing  Measurements 

The  goal  of  this  task  was  to  develop  a  measurement  capability  at  ISTEF  to  quantify  the 
atmospheric  seeing.  "Seeing^  throu^  the  atmo^here  is  degraded  by  turbulent  iluauations  in  the 
atmo^herk  refractive  index.  CREOL  successfully  demonstrated  the  use  of  an  experimental  setup 
consisting  of  a  telescope  and  a  video  camera  to  measure  the  atmospheric  seeing.  The  seeing  is 
quantified  by  the  Fried  parameter,  Using  this  CREOL-developed  instrumentation,  it  was  shown 
that  at  5  /im  and  for  "good”  seeing  conditions  at  ISTEF  and  a  near  vertical  path,  r„  was  on  the  order 
of  14  cm. 

13  Task  HL  Atmospheric  Induced  Enhanced  Backscatter  Characterization 

The  goal  of  this  task  was  to  characterize  enhanced  backscatter  (BBS),  which  occurs  to  a 
laser  beam  when  it  passes  through  the  same  turbulent  atmosphere  twice.  The  double  passage  occurs 
in  mono-static  radar  systems.  The  phenomenon  was  simulated  on  a  computer,  and  numerical 
experiments  conducted.  Laboratory  experiments  were  conducted  on  a  laser  beam  passing  twice 
through  a  random  phase  screen.  The  EBS  signal  was  found  with  simple  averaging  of  a  TV  image. 
Joint  work  was  conducted  with  the  British  MoD  Royal  Signals  and  Radar  Establishment  on  an 
outdoor  1000  meter  range  experiment.  The  averaging  of  TV  images  clearly  shows  diffraction 
limited  information  as  to  the  target  location. 

1.4  Task  rv.  IR  Focal  Plane  Array  Testing 

CREOL’s  Focal  Plane  Array  testing  program  consists  of  the  measurement  of  various  specific 
detectors,  the  development  of  various  instrumentation  devices  for  the  measurement  of  detector 
modulation  transfer  functions  (MTFs),  the  development  of  an  IR  scene  projector  system  as  well  as 
the  development  of  data  processing  and  computing  facilities. 

A  new  technique,  utilizing  a  ^>eckle  generator,  for  the  measurement  of  the  detector 
modulation  transfer  function  (MTF)  has  been  developed.  For  the  first  time  MTF  measurements 
of  hexagonal  arrays,  utilized  in  super-centroiding  applications,  have  been  reported. 


2.  ONE  MICRON  LASER  RADAR 


Philip  Gatt,  Ph.D. 

Laser  Radar  Laboratory /CREOL 

2.1  Background 

The  goal  of  this  task  was  to  provide  ISTEF  with  a  one  micron  ladar  capability  and  for  plume 
backscatter  measurements  and  for  evaluation  and  development  of  innovative  ladar  systems  and 
components.  A  coherent  one  micron  ladar  will  provide  1ST  with  a  target  discrimination  and 
tracking  capability,  providing  an  increased  (10  times)  resolution  over  that  of  the  more  mature  ten 
micron  coherent  ladar  technology.  Significant  scientific  innovations  that  relate  specifically  to  shorter 
wavelength  ladars  should  be  investigated.  These  innovations  include  the  use  of  tracking  local 
oscillator,  required  for  the  larger  anticipated  doppler  shifts,  a  fiber  mixer/coupler  to  improve  system 
ruggedness,  and  an  array  receiver  to  enhance  the  signal-to-noise  ratio  (SNR)  in  the  presence  of 
atmospheric  turbulence 

The  initial  task  of  CREOL’s  one  micron  ladar  program  was  to  mechanically  harden  an 
existing  laboratory  grade  GFE  l.OJ,  Nd:YAG  coherent  laser  radar  system  for  field  operation  in  an 
on-the-mount  config’.tration.  This  laser  radar  system  was  to  be  supplied  by  another  government 
contractor.  Due  to  difficulties  associated  with  the  delivery  of  this  government  asset,  CREOL  was 
directed  to  design  and  construct  an  alternate  prototype  laser  radar  system  to  be  deployed  at  ISTEF. 
This  coherent  ladar  ^stem  would  incorporate  a  10cm  monostatic  off-axis  transceiver  and  a 
government  furnished  1.2J  20Hz  1.06/im  laser  purchased  through  Patrick  Air  Force  Base  (PAFB). 
Purchasing  delays  at  PAFB  resulted  in  a  late  delivery  (January  1991)  of  this  laser.  During  the 
PAFB  laser  procurement  cycle,  CREOL  purchased,  fabricated,  and  tested  the  major  subsystems  for 
the  1.06^  coherent  ladar  system. 

Due  to  changing  scientific  goals  CREOL  was  redirected  in  January  of  1991  to  design  and 
fabricate  a  direct-detected  one  micron  laser  radar  rather  than  a  coherent  laser  radar,  for  the 
purpose  of  conducting  plume  backscatter  measurements  and  an  autodyne  experiment,  after  which 
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the  CREOL  coherent  laser  radar  program  was  to  be  resumed.  The  plume  backscatter  and  the 
autodyne  e]q>enments  will  utilize  a  Contraves  bistatic  10cm  Coude’  feed  transmitter  and  SOcm 
pseudo-Cassegrain  receiver  telescope  and  mount  system.  This  telescope  is  scheduled  for  ISTEF 
delivery  from  San  Di^p  in  the  fall  of  1991. 

The  autodyne  e3q)eriment  is  a  "proof  of  concept,"  active  (532nm)  target  tracking  experiment. 
This  concq)t  tracks  the  target  at  the  plume/hardbody  interface  by  maintaining  a  track  on  a  "self¬ 
beat"  note  generated  by  a  velocity  difference  between  the  hardbody  of  the  target  and  its  plume.  The 
plume  backscatter  experiment  will  employ  a  direct  detected  laser  radar  to  generate  plume  cross 
section/reflectance  data  at  both  1064  and  532  nm.  These  data  will  be  collected  for  an  ensemble  of 
"targets  of  t^jportunity*  available  to  the  ISTEF  range. 

12  Coherent  Laser  Radar  Development 

During  the  PAFB  procurement  cycle  of  the  laser,  major  subsystems  for  the  coherent 
receiver  (Figure  2-1)  were  purchased,  fabricated,  and  tested.  These  subsystems  include;  1)  an  off- 
axis  eccentric  pupil  lOX  10cm  beam  expander  to  be  employed  as  the  monostatic  transceiver,  2)  a 
single-mode  diode-pumped  5mW  CW  NdrYAG  monolithic  ring  laser  with  fast  PZT  tuning  control 
and  slow  temperature  control  to  facilitate  offset  frequency  locking  to  the  transmitter,  3)  a  2GHz  p-i- 
n  InGaAs  detector  with  a  quantum  efGciency  of  70%  and  a  linear  range  greater  than  5  mW  of 
optical  power,  4)  wide  bandwidth,  low  noise,  RF  preamplifier  and  amplifiers,  and  5)  high  power  AR 
coated  optics. 

The  beam  expander  has  a  standard  aluminum  reflective  coating  with  a  protective  SiO 
overcoat.  This  coating  is  sufficient  assuming  a  long  pulse,  10^  format;  however,  for  an  8ns  pulse 
the  power  level  (1.5  GW/cm^)  exceeds  the  damage  threshold  of  aluminum.  If  this  expander  is  to 
be  used  as  a  short  pulse  transmitter,  the  secondary  mirror  must  be  recoated  with  a  high  power  HR 
coating.  Although  it  would  be  desirable  to  have  a  broad  band  coating,  none  can  be  obtained  for 
these  high  power  q>ecillcations.  High  power,  high  reflective  dielectrics  can  easily  handle  the 
damage  threshold  for  single  wavelength  operation.  However,  these  power  levels  push  the  state  of 
the  art  for  high  power  harmonic  HR  coatings.  CREOL  has  successfully  tested  a  witness  sample  of 
such  a  harmonic  HR  coating  and  will  have  both  the  secondary  and  primary  stripped  and  recoated 
with  this  dielectric. 

All  other  components  of  the  coherent  laser  radar  were  tested  and  performed  well  within  the 
manufacturer’s  specifications.  Of  particular  interest  to  a  coherent  laser  radar  is  the  upper  limit  of 
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Figure  2-1  Coherent  Laser  Radar  Block  Diagram 


COHERENT  1.06  MICRON  LADAR  SCHEMATIC 
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the  detector’s  linear  range,  «^ch  for  the  Epitax  ETX  lOOT  was  found  to  be  well  above  5  mW  of 
optical  power.  Since  the  incoherent  ladar  program  has  superseded  the  coherent  program, 
integration  and  testing  of  these  components  have  been  delayed. 

23  Continuing  Research 

CREOL  is  currently  operating  under  a  follow-on  ONR  grant.  Under  this  grant  CREOL  will 
perform  the  following  activities  given  in  order  of  priority:  1)  perform  laser  characterization  and 
modifications,  2)  continue  the  development  of  an  incoherent  laser  radar  system,  3)  conduct  initial 
ISTEF  incoherent  ladar  eiqjeriments,  4)  train  the  ISTEF  personnel,  5)  develop  and  exercise  an  in- 
house  coherent  one  micron  laser  radar  simulation  capability,  6)  convert  the  incoherent  one  micron 
ladar  to  a  coherent  one  micron  ladar,  and  7)  transfer  the  technologies  developed  in  the  coherent 
ladar  simulation  lab  into  the  ISTEF  coherent  ladar  testbed. 

1  Laser  Characterization  and  Modifications 

Laser  characterization  is  a  critical  step  in  the  design  of  a  ladar  system.  A  characterization 
of  both  the  near  and  far  field  profiles  is  critical.  The  near-field  profile  must  be  carefully  analNzed 
to  determine  the  peak  intensity  levels  of  localized  hot  spots.  The  far-field  profile  is  critical,  as  this 
parameter  determines  the  beam  divergence,  energy  transfer  efficiency,  and  beam  shape  as  needed 
to  invert  reflectivity  data. 

It  has  been  determined  that  the  laser  (a  Continuum  NY82-20)  suffers  from  excessive 
divergence  in  the  green  and  residual  IR  beams.  Even  though  the  manufacturer  did  not  spec  the 
beam  quality  output  of  the  doubler,  from  discussions  with  a  the  laser  designer,  CREOL  has 
determined  that  the  system  is  behaving  abnormally.  CREOL  is  currently  investigating  this  problem 
and  making  the  necessary  modifications  to  correct  the  system. 


2.3.2  One  Micron  Incoherent  Ladar  Development 

CREOL  is  currently  (July  1991)  completing  the  design  of  the  transmitter  and  receiver  of  a 
field  hardened  one  micron  incoherent  ladar.  The  details  of  this  design  will  be  forthcoming  in  a  final 
report  to  be  submitted  in  January  of  1992.  Remaining  tasks  include:  final  system  assembly,  detector 
(APD)  testing  and  calibration,  and  range  gate  electronics  design. 


CREOL  will  conduct  initial  ladar  experiments  at  ISTEF  to  verify  system  operation  and 
calibration.  These  e]q>eriments  will  include,  in  chronological  order:  1)  tests  at  the  ISTEF  1  km 
range  against  stationary  targets,  2)  low  altitude  satellite  tests  (assuming  approval  is  obtained),  and 
3)  tests  against  targets  of  opportunity  including  both  plume  and  hard-bodies. 

2^i-ISIEF  Pgnonngl  TraininB 

CREOL  will  ensure  that  the  appropriate  ISTEF  (site  contractor)  personnel  are  briefed  and 
trained  for  system  operation.  This  training  will  occur  simultaneously  with  the  initial  tests  mentioned 
above.  After  system  handoS,  test,  and  training,  CREOL  will  be  on  call  for  further  ISTEF  support. 

2.3.5  One  Micron  Ladar  Simulation  Laboratory 

CREOL  is  beginning  the  development  of  a  one  micron  ladar  simulation  laboratory.  This 
laboratory  will  include  both  a  computer  simulation  capability  as  well  as  an  analog  simulation 
laboratory.  Areas  of  research  interest  include;  1)  innovative  receiver  structures  and  signal 
processing  to  improve  heterodyne  efGciency  in  the  presence  of  atmospheric  turbulence,  2) 
innovative  techniques  to  harden  a  coherent  transceiver,  3)  measurement  and  prediction  of  ISTEF 
atmospheric  extinction  conditions  and  calibration  techniques,  4)  laser  radar  cross-section  analysis 
and  measurement,  and  5)  bidirectional  laser  reflectance  data. 

2.3.6  Coherent  One  Micron  Ladar  Development 

After  the  successful  completion  of  the  plume  backscatter  and  autodyne  experiments,  the 
Nd:YAG  ladar  is  scheduled  to  be  reconfigured  into  a  coherent  laser  radar  system.  This  ladar  will 
be  a  single  frequency,  1064  nm,  employing  a  4"  monostatic  transceiver  and  operated  in  an  "on-the- 
mounf  configuration  on  the  ISTEF  Contrave.*:  Kineto  mount.  In  the  initial  configuration,  this  ladar 
will  utilize  a  single  InGaAs  p-i-n  detector  and  a  variable  offset  local  oscillator,  with  a  conventional 
beam  splitter  for  local  oscillator  mixing. 

2.3.7  Technology  Transfer 


CREOL  anticipates  modifying  the  initial  coherent  laser  radar  configuration,  described  above, 
to  employ  an  array  receiver  and  fiber  optic  coupling  technology  developed  in  the  1.06/im  Ladar 
Simulation  Laboratory  at  CREOL.  This  technology  will  allow  larger  aperture  receivers  to  be 
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employed  in  the  coherent  ladar  program,  while  maintaining  heterodyne  efficiency  and  increasing 
the  heterodyne  carrier  to  noise  ratio.  Current  plans  are  to  employ  the  4*/20"  optic  incorporated 
in  the  Coude’  feed  mount,  which  will  be  used  in  the  incoherent  Nd;YAG  e:q>eriments. 


23.8  Program  Schedule 

Below  is  an  estimated  time  schedule  for  the  continued  research  efforts  to  be  conducted  at 
CREOL. 


CREOL  Nd-YAv-  ogram  Time  Schedule 
1991  1992 

3  6  9  12  3  6 

Incoherent  System  Design  •  •  • 

Incoherent  Experiments  •  •  » 

ISTEF  Personnel  Training  * 

Simulation  Lab  «  »  »  • 

Coherent  Ladar  Development  «  •  » 

Coherent  Experiments 

Coherent  Technology  Transfer  * 


3.  ATMOSPHERIC  SEEING  MEASUREMENTS 

Arthur  R.  Weeks,  Ph.D. 

Image  Processing  Laboratory/Computer  Engineering 

3.1  Background 

The  purpose  of  this  task  was  for  CREOL  to  provide  ISTEF  with  a  measurement  technique 
which  ’atmo^heric  seeing"  could  be  quantified.  "Atmospherir  seeing"  is  defined  as  the 
degradation  of  an  image  by  the  atmosphere  due  to  atmospheric  turbulence.  Since  all  of  the  images 
completed  and  proposed  at  ISTEF  require  long  range  imaging  through  the  atmosphere,  it  is  of 
primary  importance  that  the  effects  of  the  atmosphere  be  quantitatively  defined. 

It  is  important  to  realize  that  the  atmosphere  degrades  an  image  in  two  ways.  The  first, 
known  as  scintillation,  is  the  result  of  turbulence  in  the  upper  atmosphere  (>3  km)  and  causes  the 
amplitude  of  an  image  to  fluctuate  with  time.  The  second  phenomenon,  known  as  seeing,  tends  to 
degrade  an  image  by  blurring  which  results  from  turbulence  just  above  the  telescope.  The  type  of 
imaging  required  at  ISTEF  is  typically  through  a  vertical  path  in  the  atmosphere  in  which  the  effects 
of  image  degradation  are  considered  the  most  important.  The  reason  that  seeing  is  the 
predominant  degradation  effect  is  that  most  of  the  atmosphere  is  contained  within  the  first  3  km 
distance  above  the  earth.  Another  important  fact  is  that  since  scintillation  effects  are  due  to 
atmospheric  turbulence  above  3  km  while  seeing  degradation  is  due  to  atmospheric  turbulence 
below  3  km,  there  is  very  little  correlation  between  the  two  degradations. 

Seeing  degradation  will  affect  an  image  in  two  ways:  first,  it  will  tend  to  make  an  object 
within  an  image  move  randomly  about  a  point,  and  second  it  will  result  in  blurring  of  the  object. 
How  these  two  affect  an  image  depends  on  the  integration  time  of  the  imaging  system.  For  long 
time  exposures,  image  motion  integrates  to  a  time-averaged  image  that  is  blurred.  Blurring  is  due 
to  many  correlation  cells  (seeing  disks)  that  appear  at  the  entrance  aperture  of  a  large  telescope 
system. 
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Each  seeing  disk  effectively  images  into  the  focal  plane  of  a  telescope  an  image  that  is 
slightly  di^laced  from  other  seeing  disks.  The  sum  of  all  of  these  seeing  disks  is  the  observed 
image,  which  is  blurred. 

In  Figure  3-1,  the  aperture  of  the  telescope  is  such  that  the  three  seeing  disks  image  an 
object  into  the  focal  plane  of  the  telescope  at  different  locations.  The  net  effect  is  a  blurring  of  the 
object.  If  the  telescope  aperture  is  reduced  in  size  so  that  only  one  seeing  disk  images  through  the 
telescope,  there  will  be  only  one  image  that  will  be  moving  randomly  about  a  fixed  central  point. 
The  random  motion  is  due  to  movement  of  the  seeing  disk  with  time.  This  is  why  when  very  short 
exposures  are  taken  of  images  observed  through  the  atmosphere,  near  diffraction  limited  results  can 
be  obtained.  A  short  timed  cjqjosure  freezes  the  image  motion  so  that  effectively  only  one  image 
is  obtained.  On  the  other  hand,  long  integration  times  integrate  this  image  motion  into  a  blurred 
unresolved  image. 

32  TTieory 

Figure  3-1  shows  that  the  size  of  the  seeing  disk  determines  the  effective  aperture  size 
needed  for  a  telescope  system.  Hence  by  measuring  the  size  of  the  seeing  disk  at  various 
observation  sites,  relative  seeing  performance  can  be  made  (Ref.  2,  4,  5),  The  parameter  used  to 
measure  the  seeing  disk  size  is  known  as  the  Fried  parameter,  r„  (Ref.  1). 

To  measure  the  Fried  parameter,  a  telescope  system  images  an  unresolvable  star,  and  the 
effeaive  size  of  star  is  determined,  assuring  the  model  for  the  star  is  a  point  source. 

The  maximum  theoretical  angular  resolution  0  of  a  telescope  is  given  by  9  =  1.2  X/D,  where 
X  is  the  wavelength  of  light  assumed  to  be  in  the  middle  of  the  visible  band  and  D  is  the  diameter 
of  the  telescope.  For  typical  values,  D  =  11  inch  and  X=0.5  ^m,  ©=  1.78  firad.  This  value  is  about 
a  factor  of  5  smaller  than  the  minimum  angular  spread  of  the  star’s  image  by  the  atmosphere.  A 
typical  spread  of  the  atmosphere  is  on  the  order  of  10  /itad,  reducing  the  effective  diameter  of  the 
telescope  to  approximately  2  inches.  This  would  tend  to  indicate  that  a  20  inch  telescope  yields  the 
same  resolution  as  a  2  inch  telescope. 

This  effective  diameter  of  the  telescope  as  mentioned  above  can  be  determined  from  Fried’s 
parameter,  r^,  which  is  a  measure  of  the  spatial  coherence  of  the  EM  wave  due  to  turbulence.  The 
intensity  pattern  at  the  image  plane  of  the  telescope  due  to  seeing  effects  is  very  similar  to  a 
Gaussian  shape.  To  find  r^,  the  Modulation  Transfer  Function  (MTF)  of  the  atmosphere  must  be 
found.  The  MTF  of  the  image  is  defined  as  the  Fourier  transform  of  the  image  of  a  point  source, 
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Figure  3-1  Typical  Telescope  System  Depicting  Multiple  Seeing  Disks 


at  the  image  plane  of  the  telescope.  The  total  MTF  of  a  star  that  is  imaged  by  a  telescope  is  given 
by 


MTF,  -  (MTF.)  (MTF.)  (MTF,).  (3-1) 

Where  MTF,,  is  the  MTF  of  the  image;  MTF„,  is  the  MTF  of  the  star;  MTF„  is  the  MTF  of  the 
atmosphere;  and  MTF„  is  the  MTF  of  the  telescope.  For  a  high  quality  telescope,  its  MTF  is 
^atially  flat  (ie.  a  constant).  A  star  which  is  unresolvable  by  the  telescope,  has  a  spatially  flat  MTF 
as  well  Thus,  the  total  MTF  of  the  imaging  ^stem  depends  only  on  the  atmospheric  MTF,  which 
is  given  by 


.  (3-2) 

Where  p  -  Xf*-;  X,  is  the  wavelength  of  light;  f,  is  the  focal  length  of  telescope;  and  y,  is  the  spatial 
frequency  in  lines/nm. 

Since  the  MTF  of  the  image  is  nearly  Gaussian  in  shape  (i.e.,  5/3  instead  of  2),  the  image 
intensity  distribution  will  also  look  very  similar  to  a  Gaussian.  The  larger  the  measured  value  of 
r„,  the  better  the  seeing  conditions  are. 

The  determination  of  r.  can  be  used  to  find  optimum  locations  or  sites  to  locate 
observatories.  Experimental  measurements  can  be  made  with  known  diameter  stars  to  determine 
the  value  of  for  various  locations.  In  addition,  the  effects  of  humidity,  temperature,  and  such 
things  as  the  height  of  the  telescope  above  the  ground  on  r„  can  be  determined. 

To  measure  directly,  the  MTF  of  the  image  is  obtained  and  the  e'^  *^  point  is  determined; 
unfortunately  this  measurement  point  is  located  within  the  low  signal-to-noise  portion  of  the  MTF 
curve,  far  out  in  the  tails  of  the  measured  intensity  distribution.  Tutomirski  and  Yura  proposed 
using  (Ref.  3) 


MTF.(p)  - 


(3-3) 


where  =  2.\p„  to  find  the  e  *  point.  This  leads  to  more  accurate  results,  since  the  measurement 
point  is  now  located  at  the  e  ’  point  instead  of  the  e’^'**  point. 
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3  J  ExperiiNat 

Figure  3-2  shows  the  experimenul  setup  used  at  ISTEF  to  measure  r,.  The  experimental 
sett^  oonUined  a  Questar  telescope,  a  Sony  CCD  camera,  the  Contraves  tradcer  system  located  at 
ISTEF,  and  an  IBM  PC  done  with  a  DaU  Translation  imaging  system.  The  Questar  telescope  is 
a  Cassagranian  telescope  with  a  12  inch  aperture  and  180  inch  focal  length.  Attached  to  the 
telescope  is  the  Sony  CCD  camera  with  a  horizontal  and  vertical  pixel  resolution,  respectively,  of 
23  by  13.4  The  video  output  of  the  camera  is  next  digitized  by  the  IBM  PC  frame  grabber 
system.  This  system  provides  the  capability  of  averaging  at  a  rate  of  8  frames/second  up  to  256 
frames. 

The  tracker  is  used  to  bring  the  image  of  the  desired  star  into  the  FOV  of  the 
tdescope/camera  imaging  system.  Once  in  the  FOV  of  the  system,  the  tracker  is  locked  into  the 
star  and  tracks  it.  When  data  are  to  be  taken,  the  star  is  moved  to  the  upper  left  corner  of  the 
viewing  ^tem,  by  applying  an  offset  to  the  tracker  system.  The  tracker  is  then  turned  off,  and  the 
star  is  allowed  to  drift  freely  through  the  FOV  of  the  telescope.  During  this  time  the  IBM  PC 
system  acquires  a  set  of  integrated  frames.  Turning  off  the  tracking  system  guarantees  that  the 
effects  of  tracker  jitter  will  be  completely  eliminated  from  influencing  the  experimental  data. 

Since  the  video  system/framc  grabber  is  an  RS-170  8  bit  system,  the  total  integration  time 
in  which  data  were  acquired  was  set  to  20  frames.  This  reduces  the  effect  of  quantization  noise 
when  computing  r„  from  the  acquired  data. 

With  this  setup,  the  procedure  required  to  obtain  the  desired  image  data  is  (see  Appendix 
A  for  more  details); 

1.  Set  up  equipment  as  shown  in  Figure  3-2. 

2.  Make  sure  that  all  video  lines  are  property  terminated.  This  reduces  video  smearing. 

3.  Use  the  IBM  PC  with  its  frame  grabber  boards  to  acquire  20  frames  of  data  per  run. 

4.  Save  these  data  on  floppy  disk  for  further  analysis. 

Performing  these  steps  will  produce  a  digitized  image  that  contains  an  integration  of  the  motion  of 
the  star  as  it  moves  through  the  FOV  of  the  telescope. 

3.4  Results 

The  combined  camera  and  telescope  system  will  image  the  star  onto  the  CCD  with  a 
diameter  D  with  the  diameter  being  proportional  to  the  focal  length  of  the  telescope.  The  total 
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Figure  3-2  CREOL  Designed  ISTBF  Seeing  Measurement  Apparatus 
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diameter  of  the  image  on  the  CXD  detector  wilJ  be 

D  ~fe  (3-4) 

where  6  is  the  angular  width  of  the  star’s  line  spread  function. 

For  6  of  about  2  arcseconds  of  seeing  and  /  «  180  inches,  the  size  of  the  star  at  the  image 
plane  will  be  45  This  is  about  2-3  pixels  of  width. 

Since  a  simple  esqperimental  setup  is  required,  the  use  of  no  tracker  reduces  the 
complications  caused  by  vibrations  from  the  telescope’s  motion,  which  blurs  the  time  averaged  line 
^read  function  even  more.  The  resulting  image  on  the  CCD  camera  will  be  a  line  produced  by  the 
star  as  it  traverses  the  telescope’s  FOV  and  averaged  by  the  CCD  camera.  The  transverse  width 
of  this  line  wiU  be  the  one-dimensional  line  spread  function. 

The  averaged  line  spread  function  is  obtained  by  averaging  several  perpendicular  slices.  The 
total  averaging  time  required  to  produce  an  averaged  line  spread  function  is  on  the  order  of  tenths 
of  second.  To  obtain  the  one-dimensional  line  spread  function,  several  transverse  slices  are  made 
through  the  horizontal  star  streak,  and  then  these  sUces  are  averaged.  Before  actually  averaging 
the  slices,  to  remove  any  low-frequency  vibrations  that  might  be  due  to  wind  vibrations  on  the 
telescope,  the  centroid  for  each  slice  is  computed.  Next,  each  shce  is  aligned  to  its  centroid  and 
averaged  with  all  the  other  slices  to  yield  the  time-averaged  one-dimensional  line  spread  function. 

To  obtain  the  MTF  of  the  atmosphere,  the  Fourier  transform  of  the  one-dimensional  line 
spread  function  is  computed  and  the  location  of  the  e  ’  point  is  determined.  Then,  with  the  help 
of  Eqn.  3-3,  can  be  computed.  It  should  be  mentioned  that,  if  the  aperture  of  the  telescope  is 
on  the  order  of  r,,  the  MTF  of  the  atmosphere  must  be  corrected  to  include  the  effects  of  the 
telescopte.  This  is  done  by  dividing  the  measured  MTF  by  the  MTF  of  the  telescope.  For  the 
experimental  setup  utilized  at  ISTEF,  the  telescope  MTF  was  considered  constant,  since  on  a  night 
of  very  good  seeing  was  typically  smaller  than  150  mm,  which  is  smaller  than  the  diameter  of  the 
telescope,  r,  can  then  be  determined  by  simply  measuring  the  e'^  points  of  the  line  spread  function 
and  computing  r„  directly.  This  step  eliminates  the  need  to  use  a  discrete  fourier  transform  to 
compute  r,  making  the  data  analysis  simpler  (Refs.  2,  4,  5). 

On  one  particular  night,  the  seeing  was  extremely  good.  The  quality  of  the  images  of  the 
stars  produced  exceptional  good  quality  even  to  the  naked  eye.  The  star  Alpha  Lyrae  was  observed, 
which  was  about  15  degrees  from  the  zenith. 
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A  curve  fit  of  this  line  ^read  function  to  a  Gaussian  curve  is 


L(x)  -  140^-“-“*)*'“, 


(3-5) 


from  tKfrkh  we  obtain 


a-ZlTd, 


(3-6) 


where  d«23  /on.  Taking  the  Fourier  transform  of  Eqn.  3-5  and  equating  it  to  Eqn.  3-3  yields 


P,-0.67 


hi 

d 


(3-7) 


or 


r  -1.41 

C 


d  ■ 


(3-8) 


Assuming  X  is  in  the  middle  of  the  visible  spectrum,  X=0.5  ftTn,f=  180  inches,  and  d=23  nm,  is 
given  by  140. 1  mm. 

To  compute  the  effective  angular  resolution,  the  diffraction  equation  (©  =  1.22  X/d)  is  used 
for  the  case  d-r^  which  yields  0=1  arcsec,  considered  to  be  very  good  seeing. 


3.5  Conciosion 

Data  were  taken  to  show  that,  using  a  simple  experimental  technique  of  a  telescope  and  a 
video  camera,  seeing  measurements  can  be  easily  made.  Using  a  portable  tracker,  which  is  now 
available  at  ISTEF,  this  equipment  setup  can  be  easily  made  portable,  so  that  remote  site  testing 
can  be  done. 

3.6  Appendix  3-A 

Equipment  Procedures: 

1.  Mount  the  Questar  telescope  to  the  Contraves  tracker  mount. 

2.  Attach  the  Sony  CCD  camera  to  the  Questar  telescope. 

3.  Attach  the  Motorola  computer  to  the  Contraves  tracker. 

Install  the  data  translation  boards  in  the  IBM  PC  clone. 


4. 
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5.  Connect  the  video  cables  from  the  Sony  CCD  camera  to  the  IBM  PC  computer  and  to  a 
video  monitor  for  viewing. 

6.  Terminate  the  video  cables  correcdy. 

7.  Run  the  IRIS  program  in  the  PC  to  display  and  acquire  the  video  images. 

8.  Acquire  the  desired  star  with  the  tracker. 

9.  Offset  the  star  in  the  video  display  to  the  top  of  the  screen. 

10.  Turn  the  tracker  off. 

11.  Acquire  the  image  data  for  X  frames. 

12.  Save  the  image  data  to  disk. 

13.  Rq>eat  step  8  until  done. 
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4.  AmOSPHERIC  INDUCED 
ENHANCED  BACKSCATTER  CHARACTERIZATION 


Ron  Phillips,  Ph.D. 

Wave  Propagation  Laboratory/CREOL 

4.1  Background 

The  goal  of  this  task  was  to  investigate  the  scattering  of  an  optical  wave  due  to  the  double 
passage  through  atmospheric  turbulence  known  as  the  enhanced  backscatter  effect.  TTie 
investigation  was  concerned  with  developing  an  understanding  of  the  physical  phenomenon  via  both 
numerical  and  laboratory  experiments.  Possible  applications  for  tracking  targets  from  the  ISTEF 
facility  were  to  be  considered. 

To  extract  diffi'action  limited  information  an  adaptive  optical  system  operates  as  a  closed 
loop  system,  whereby  an  error  signal  of  a  physical  observable  is  minimized.  If  an  optical  system 
were  to  operate  in  an  open  loop  mode  and  still  extract  diffraction  limited  information,  a  most 
sophisticated  understanding  of  the  propagation  process  would  be  required.  A  closed  loop  system 
simply  minimizes  an  error  by  physically  changing  the  optical  system  such  as  by  a  deformable  mirror. 
An  open  loop  system  would  be  a  much  simpler  hardware  system  but  would  require  signal 
processing,  thereby  putting  the  burden  of  the  information  extraction  on  computation  and  NOT  on 
hardware.  Usually  this  is  a  less  expensive  solution,  since  electronic  hardware  is  less  expensive  than 
optical  hardware.  The  open  loop  approach  is  model  dependent  and  requires  a  simple  yet 
sophisticated  mathematical  model  if  the  information  is  to  be  extracted  in  real  time.  Hence,  open 
loop  op>eration  of  an  optical  system  that  can  extract  diffraaion  limited  information  requires  a  deep 
understanding  of  the  propagation  physics. 

There  are  two  basic  architectures  for  laser  beam  systems  for  tracking  and  imaging  of  targets. 
The  bi-static  system  is  an  architecture  with  the  laser  transmitter  aperture  separate  from  the  receiver 
aperture.  TTie  laser  beam  path  can  be  assumed  to  be  through  two  independent  turbulent 
propagation  paths.  The  second  architecture  is  a  mono-static  system,  in  which  the  laser  transmitter 
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and  receiver  share  the  same  output  aperture.  Hence  the  transmitter  and  receiver  have  coaxial  fields 
of  view.  The  propagation  path  of  the  laser  beam  then  passes  twice  through  the  same  turbulence 
(assuming  slow  evolution  of  the  turbulence).  The  scattering  of  the  returned  wave  is  then  correlated 
with  the  tran»nitted  wave.  This  correlation  permits  a  reciprocity  of  the  passage  of  a  wave  through 
any  scattering  path.  In  other  words,  two  rays  can  pass  through  any  scattering  path  in  opposite 
directions  and  emerge  in  phase  with  one  another.  It  is  the  presence  of  correlation  that  may  be 
exploited  to  extract  diffraction  limited  information  from  the  returned  beam. 

The  physical  genesis  of  the  phenomenon  of  enhanced  backscatter  (EBS)  is  reciprocity  of 
electromagnetic  wave  propagation.  If  a  transmitted  ray  is  reflected  from  the  target  back  to  the 
receiver,  another  ray  can  take  the  same  path  but  in  the  reverse  direction.  The  two  rays  will  arrive 
at  the  receiver  in  phase  with  one  another  because  they  both  traversed  the  same  inhomogeneities, 
albeit  in  reverse  order. 

Depending  upon  the  specifics  of  the  wave,  i.e.,  plane,  spherical,  or  beam  wave,  the  double 
pass  scattering  phenomenon  has  the  following  characteristics: 

1.  An  increase  in  the  fluctuations  of  the  intensity  at  the  transmitter  aperture  plane; 

2.  An  increase  in  the  fluctuations  of  the  phase  of  the  wave; 

3.  A  decrease  in  the  transverse  spatial  phase  of  the  wave; 

4.  A  retention  of  diffraction  limited  information  about  the  reflecting  target. 

42  Analysis  of  Enhanced  Backscatter 

The  essence  of  EBS  can  be  understood  by  considering  a  uniform  plane  wave  passing  through 
a  random  phase  screen  and  being  reflected  back  onto  itself.  Here  we  use  a  simple  phase  screen  as 
a  model  for  the  turbulent  atmosphere  and  a  ray  model  for  the  wave.  For  a  uniform  plane  wave 
incident  on  the  screen,  the  rays  are  all  parallel  and  propagating  in  the  direction  of  incidence.  The 
medium  disperses  the  rays  over  a  variety  of  directions.  The  rays  then  propagate  onto  the  target  and 
are  reflected  back  through  the  same  screen.  The  rays  emerging  from  the  screen  along  reciprocal 
paths  will  emerge  in  phase  in  the  direction  of  incidence  of  the  incoming  wave.  Although  the 
absolute  phase  of  the  scattered  wave  will  fluctuate,  in  the  direction  of  the  incidence  the  scattered 
waves  will  be  in  phase  because  the  rays  traverse  the  same  refractive  inhomogeneities  but  in  opposite 
directions.  In  Figure  4-1,  the  rays  Ej  and  are  in  phase  and  Ej  and  Ej  are  in  phase,  but  Ej  and 
Ej  are  not  in  phase,  nor  are  Ej  and  E4,  because  these  rays  do  not  have  reciprocal  paths.  Suppose 
we  now  place  a  coUection  lens  behind  the  aperture  and  an  array  detector  behind  the  lens.  If  the 
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array  detector  is  in  the  focal  plane  of  the  lens,  it  will  be  in  the  far  field  of  the  aperture.  The 
detected  double  pass  coherent  wave  will  appear  as  a  sharp  intensity  peak  in  the  focal  plane.  The 
location  of  this  peak  relative  to  the  entire  received  signal  will  give  the  location  of  the  target  to  the 
resolution  of  the  optical  system. 

4  J  Ejqferimental  Results 

Figure  4-2.  is  a  series  of  photographs  of  the  intensity  fluctuations  of  a  laser  beam  after 
passing  through  a  random  phase  screen  generated  by  warm  air  turbulence.  Figure  4-2a  is  just 
beyond  the  screen,  Figure  4-2b  is  at  the  "focus"  of  the  screen,  and  Figure  4-2c  is  in  the  far  field  of 
the  screen.  Figure  4-3  shows  a  TV  image  of  a  plane  wave  after  passing  twice  through  a  random 
phase  screen  in  a  laboratory  experiment.  Figure  4-3a  shows  the  characteristic  speckle.  Figure  4-3b 
is  the  TV  image  of  the  multi-frame  averaging  automatically  performed  by  the  TV  camera  by  rapidly 
moving  the  random  phase  screen  faster  than  the  camera  could  time  resolve  the  moving  speckle 
pattern.  The  bright  spot  is  the  enhanced  backscatter  due  to  the  double  passage. 

Figure  4-4  is  a  photograph  of  the  laser  intensity  fluctuations  after  propagating  a  distance  of 
1000  meters  on  an  outdoor  range.  The  aperture  of  the  collecting  lens  for  the  camera  was  smaller 
than  the  beam  diameter. 

Joint  experiments  were  performed  with  the  Royal  Signals  and  Radar  Establishment  of  the 
British  Ministry  of  Defense  in  Malvern  England.  EBS  experiments  have  been  performed  on  a  1,000 
meter  outdoor  range  (see  Figure  4-5a).  A  HeNe  laser  beam  was  pointed  down  range  to  a 
retro-reflector.  The  returned  wave  was  detected  by  a  CCD  TV  camera  located  in  the  focal  plane 
cf  a  collection  lens.  The  TV  images  were  recorded  on  a  VCR.  The  TV  recorded  images  were 
played  into  a  frame  grabber  and  128  frames  averaged  together.  The  result  is  shown  in  Figure  4-5b. 
The  angular  resolution  of  the  position  of  the  target  (retro-reflector)  due  to  turbulence  is 
qualitatively  defined  by  the  smeared  signal.  The  EBS  signal  clearly  shows  the  location  of  the  target 
to  within  the  diffraction  limit  of  the  collection  optics.  The  picture  shows  an  approximate  factor  of 

5  increase  in  angular  resolution  over  the  turbulence  reduced  resolution. 

4.4  Conclusions 

The  results  of  this  study  show  that  the  enhanced  backscatter  of  the  double  passage  of  a  wave 
through  the  atmosphere  can  be  used  to  obtain  diffraction  limited  information  about  a  target.  This 
was  shown  by  processing  data  from  both  laboratory  and  outdoor  experiments. 
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Figure  4-2.  The  top  photo  is  the  intensity  fluctuation  of  a  laser  beam  just  beyond  the  random 
phase  screen.  The  middle  photo  is  the  beam  fluctuations  in  the  "focus"  region  of  the  screen. 
The  lower  photo  in  the  far-field  of  the  screen  showing  the  Gaussian  like  fluctuations  of  the 
beam. 


Fieure  4-3.  This  is  a  TV  image  of  a  beam  after  passing  tv,’ice  through  a  phase  screen.  Figure 
4-3a  IS  the  speckle  without  processing  the  signal.  Figure  4-3b  shows  the  EBS  signal  after 
averaging  many  TV  frames. 
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EBS  signal  gives  D.L. 
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location  angle 
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Uncertainty  in  target 
location  due  to 
turbulence 


Figure  4-5a.  A  schematic  of  a  mono-static  transmitier/receiver.  The  EBS  signal  can  be 
extracted  from  the  focal  plane  of  receiving  collection  lens  by  processing  the  detected  signal 
from  a  CCD  array  in  the  focal  plane. 
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A  timulation  wascnnstrurled.and  numerical  ezperimentaivere  conducted,  of  optical-wave  dnublrpaaa  propaeatinn 
thrutifh  a  phaae  acreen.  The  phaae  screen  had  a  Gauiaian  phase  fluctuation  with  a  Koimof orov  power-law  tpatial 
spectrum.  The  results  are  compared  with  a  theoretical  model  developed  for  a  aimilar  scattering  process  but  for  a 
phase  screen  with  a  Gaussian  spatial  spectrum. 


1.  INTRODUCTION 

There  has  been  a  vested  interest  in  the  study  of  waves  re¬ 
flected  or  scattered  backward  through  a  random  inhomogen¬ 
eous  medium.  Both  theoretical  and  experimental  work  has 
been  conducted  in  this  area  moat  recently  by  Jakeman,' 
Jakeman  et  al.,’  Tapster  et  ol.,^  Tatarski/  and  Eiapipov  and 
Zosemov.^  The  interest  comes  about  because  of  the  areas  of 
research  that  can  be  developed  with  this  phenomenon. 
These  areas  include  imaging,  remote  sensing,  and  the  study 
of  sound-wave  propagation  through  the  ocean.  This  en¬ 
hancement  might  be  used  to  track  objects  in  the  atmosphere. 
Furthermore,  this  phenomenon  might  be  used  to  find  the 
optical  axis  in  a  complex  system. 

The  study  of  wave  propagation  through  random  media 
followrs  the  analysis  tnat,  as  a  wave  passes  through  statisti¬ 
cally  independent  phase  screens,  the  screeru  raonotonically 
increase  the  scattering  effects  of  the  incoming  wave.  This 
increase  does  not  occur  when  a  wave  is  reflected  back 
through  the  screen  that  it  initially  uaversed.  The  effect  of 
double  passage  on  the  wave  through  the  same  inhomogene¬ 
ities  of  the  medium  is  an  enhancement  of  the  average  inten¬ 
sity  of  the  backward-reflected  wave.  This  effect  was  ob¬ 
served  in  a  review  paper  by  Kravtsov  and  Saichev.* 

A  numerical  experiment  was  conducted  on  the  phenome¬ 
non  of  enhanced  backscatter.  This  research  was  based  on  a 
simulation,  which  was  developed  by  Filice,*  of  interplane¬ 
tary  scintillation  of  radio  waves.  This  simulation  begins 
with  a  monochromatic  plane  wave  passing  through  a  random 
phase  screen.  The  wave  is  propagated  some  given  distance 
past  the  focal  plane  of  the  screen  and  reflected  back  toward 
the  same  screen.  The  average  intensity  statistics  are  stud¬ 
ied.  Although  there  are  some  limitations  to  simulations  of 
this  type,  they  do  provide  a  useful  tool  in  the  analysis  of 
various  case  studies.'-^  For  instance,  the  simulation  enables 
us  to  study  different  types  of  mirrors,  which  are  used  as 
reflectors,  at  various  angles  and  contours.  This  simulation 
has  been  tested  and  compared  with  theoretical  results  pub¬ 
lished  by  Jakeman.' 

Section  2  develops  the  physical  background  needed  for  the 
analysis,  which  is  followed  in  Section  3  by  theoretical  mod¬ 
els  Section  4  describes  the  algorithm  used.  Section  5  is  a 
brief  description  of  some  of  the  results  obtained  from  the 
simulation  Finally,  in  Section  6  we  compare  these  results 
of  the  intensity  statistics  with  theoretical  outcomes  pro¬ 
duced  by  Jakeman. 


2.  BACKGROUND 

When  an  optical  wave  passes  through  a  phase  screen,  the 
emerging  wave  has  a  perturbed  spatial  phase.  As  the  wave 
continues  to  propagate  beyond  the  screen,  the  wave  begins 
to  exhibit  amplitude  fluctuations  if  the  impressed  perturba¬ 
tions  are  large  (variance  of  phase  »  1).  These  amplitude 
fluctuations  initially  increase  with  distance:  they  reach  a 
maximum  and  then  diminish  with  further  increase  in  dis¬ 
tance  and  eventually  approach  a  saturated  value.  The  dis¬ 
tance  at  which  these  peak  fluctuations  occur  is  referred  to  as 
the  focus  distance.  These  amplitude  fluctuations  are  due  to 
two  different  mechanisms,  refractive  scattering  and  diffrac¬ 
tive  scattering.  The  relative  amount  of  each  of  these  two 
mechanisms  depends  on  the  spatial  spectrum  of  the  refrac¬ 
tive  inhomogeneities  and  the  distance  that  the  wave  has 
propagated  beyond  the  screen.  Diffractive  scattering  domi¬ 
nates  the  intensity  scintillation  for  scales  of  the  spatial  inho¬ 
mogeneities  within  the  screen  of  the  order  of  the  Fresnel 
scale,  \X2,  or  smaller.  Strong  refractive  inhomogeneities 
within  the  phase  screen,  which  are  large  compared  with  the 
Fresnel  scale,  refractively  scatter  the  wave.'" 

The  most -often-used  spatial-spectrum  models  of  the  re¬ 
fractive  inhomogeneities  are  Gaussian  and  power  law.  The 
Gaussian  spectrum  has  the  characteristic  of  smooth  inhomo¬ 
geneities  with  a  single  scale.  If  the  single  scale  of  the  Gauss¬ 
ian  spectrum  is  larger  than  the  Fresnel  scale,  the  scattering 
is  due  to  refractive  scattering  as  is  the  case  close  to  the  phase 
screen.  At  long  distances  from  the  screen,  when  the  F resnel 
scale  is  large  compared  with  the  characteristic  scale  of  the 
spectrum,  the  scattering  process  becomes  dominated  by  dif¬ 
fraction.  Hence  the  scattering  processes  are  somewhat  se¬ 
quential,  with  the  refractive  process  acting  close  to  the 
screen  and  the  diffractive  process  aaing  at  long  disunces 
from  the  screen. 

The  power-law  spectrum  has  a  range  of  inhomogeneiiies. 
Usually  the  spectrum  is  assumed  to  have  a  smallest  and  a 
largest  scale,  called  the  inner  and  outer  scales,  respectively. 
Since  the  power-law  spectrum  is  an  inverse  power  la*,  the 
power  in  the  spatial  scales  decreases  as  the  scale  gets  smaller, 
i.e..  the  spatial  wave  number  increases.  Ai  every  distance 
beyond  the  phase  screen  where  the  Fresne  -calejjes  be¬ 
tween  the  inner  and  outer  scales,  i  e.,  where  ..  < 
there  are  always  some  scales  of  inhomogeneities  contribus- 
ing  to  both  the  diffractive-  and  refractive-scattering  prtxress- 
es.  The  dominant  mechanism  depends  on  the  relative 
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■upi^rtb*  ol  the  inhooiafcncity  acalM.  At  loni:  dMUnoa* 
fruiD  the  tcrMH  where  the  FreeneJ  ecalc  caoeedt  the  outer 
•gaW.  dillrecUon  dominetea.  At  dieiencak  where  the  Free- 
■el  acAte  freetly  eicAede  the  outer  eatle,  Le.,  at  »  Lo,  the 
oaotribution  to  ooch  point  in  the  ecattered  wave  oanaa  from 
■Mny  indopoodent  inhnmoponeitiaa  and  from  the  centra) 
liaait  thoorcia:  the  field  appruachaa  oomplex  Gauaaian.  If 
tha  acraan  ia  apatially  hooioganaoua,  the  field  ia  ctrcular- 
eonaplei  Gauaaian  and  becofnea  apeckle. 

A*  the  wave  propairatas  beyond  the  phaae  acreen,  the 
tranaverae  apatia)  corrclatimi  diminiahea."  If  a  aecond 
acraan  ia  pla^  in  the  path  of  the  aravc,  the  affact  that  it  haa 
on  the  wave  dependt  on  the  wave'a  tranaverae  correlation. 
When  the  tranaverae  correlation  of  the  incident  wave  ia  lar(- 
er  than  the  refractive  inhomoiteneitiet  in  the  aecond  acreen, 
the  inhoiaapeneitiaa  contribute  to  the  phaae  fluctuation*  of 
the  wave.  For  a  tranaverae  correlation  aifnificantly  amailer 
than  the  inhotnogeneitiaa,  the  additional  effect  of  the  aecond 
acreen  ia  amall. 

When  a  e/ave  paaan  through  a  phaae  acreen  and  ia  reflect¬ 
ed  back  throufh  the  came  acreen,  the  effect  of  the  aecond 
can  five  riae  to  a  aifnificantly  different  phenomenon. 
Tliia  phenomenon  ia  an  enhancement  of  the  backacattered 
wave  and  it  due  to  the  refractive  fluctuationi  of  the  acreen 
beinc  correlated  with  the  wave.  There  are  two  aeperaie 
mechaniamt  that  can  irive  an  enhancement  in  the  backkcat- 
ter:  coherent  and  incoherent  The  coherent  mechaniam  la 
the  cofiatrucuve  interference  of  portion*  of  the  reflected 
wave  after  it  hat  travemed  different  refractive  inhomopenei- 
tie*.  Thif  It  the  dominant  contributing  mechaniam  Ui  the 
backacatter  when  the  mirror  it  placed  well  lieyond  the  focua- 
in{  diatance  of  the  phaae  acreen.  The  incoherent  mecha¬ 
niam  occur*  when  a  portion  of  the  reflected  wave  pataet 
throufh  the  aame  refractive  lenalike  inhomopeneity  that  fo- 
cuaed  thit  portion  of  the  wave  onto  the  mirror.  The  en¬ 
hancement  factor  fur  the  incoherent  caute  may  be  much 
larger  than  that  fur  the  coherent  cate,  which  it  only  a  factor 
of  2  above  the  normal  tcaiierinp  level. 

The  angular  diauibolion  of  the  backacatter  enhancement 
i*  very  narrow  for  the  coherent  procea*.  The  enhancement 
width  dacraaaaa  as  the  angular  diatribution  of  raya  increase* 
at  a  result  of  the  firat  pau  through  the  phase  acreen.  The 
width  of  the  enhancement  for  the  incoherent  process  is 
greater  than  that  of  the  coherent  process.  The  width  of  this 
contribution  ia  comparable  with  the  width  of  the  specular 
fwpiUr  distribution  of  normal  scattering. 

3.  THEORETICAL  MODELS 

A  wave  pattn  through  the  scattering  medium  and  it  reflect- 
•d  back  onto  itself.  The  backacatured  wave  emerges  from 
the  medium  into  a  eon*  of  directions.  The  emerging  wave* 
add  together  preferentially  in  the  revert*  direction  of  their 
direction  of  incidence.  For  an  esplanation  of  the  phvtical 
acattering  process,  a  ray  model  it  most  useful. 

For  a  uniform  plant  wave  incident  upon  the  screen,  the 
rays  are  all  parallel  and  propagate  in  the  direction  of  inci¬ 
dence.  Within  the  medium,  the  rays  are  dispersed  over  a 
variety  of  random  directioru.  After  pasting  through  the 
phase  screen,  the  scattered  rays  propagate  onto  the  reflect¬ 
ing  surface  and  arc  reflected  back  through  the  same  screen. 
We  assume  that,  during  the  lime  period  during  which  the 
rays  are  propagating  to  the  reflector  and  back,  the  refractive 


inhotnogencitia*  within  the  screen  do  not  change.  W'c  as- 
aume  that  the  refractive  inhomogeneitia*  change  slowly  ao 
that  aubaaquant  waves  and  their  rays  will  expermnee  other 
laaliiationa  of  the  random  inbomoganaitiaa. 

Cooaidar  aame  particular  ray  Ei  incident  upon  the  phaae 
acreen.  Now  auppoee  that  £i  emerges  from  the  aaoood  pass 
through  the  acreen  at  aome  new  location  but  pointed  in  the 
reverse  direction  of  the  incident  wave  as  shown  in  Fig.  1.  At 
the  csit  location  there  it  a  second  ray  E,  enuring  the  phaae 
screen.  This  aecond  incident  ray  will  traverse  the  same  g>ath 
at  the  first  ray  but  in  the  oppusiu  direction.  Ray  Et  will 
then  emerge  from  the  second  pass  through  the  screen  at  the 
location  where  £|  enured  the  screen  but  pointed  in  the 
uppoeiu  direction  U  £|.  The  two  ravt  eschange  locations 
and  exit  parallel  to  each  other,  because  these  two  rays 
traversed  the  same  path  but  in  oppoaiu  directions,  they  exit 
with  DO  phaae  difference  due  to  the  acaturing  process.  In 
aubsequent  raalixations  this  pairing  of  the  returning  rays 
will  occur  whenever  any  ray  can  be  scatured  back  into  the 
direction  of  the  incident  wave. 

If  *ve  coruider  rays  emerging  from  the  second  pas*  through 
the  phaae  screen  into  directions  other  than  the  incident 
direction,  there  ia  no  reciprocal  path  poasible.  Hence  there 
will  be  a  random  phase  difference  between  any  two  rays  in 
subsequent  realixation  of  the  screen  because  each  ray  will 
always  pas*  through  a  different  path.  Therefore  all  ray*  will 
be  incoherent  in  all  other  directions  but  the  direction  of 
incidence. 

Mathematically  we  can  represent  this  phx'sical  ray  model 
by  the  following  analysis.  Let  the  total  field  lie  repreaenud 
by  four  component  ra.va  such  that  the  lout  field  is 

The  inunsity  of  the  field  afur  the  double  paa*  it 

;  -  !£!■  -  1£,!-  +  1£,F  +  l£:.l-  +  l£.F  +  2£,£,  cos  f,. 

+  2£|£:,  cos  e,,  +  2£,£,  cos  f,,  *  2£.£.i  cot  e^,.  (2) 

Suppose  w*  have  a  uniform  wave  ao  that  A.  “  l£|F  “  l£}*’  “ 
lEjl-  •  l£4l-  and  therefore 

/  »  4/i,  ■+  2/i  coa  f  j-  2/i,  cos  +  27i,  cos  f  ,3  +  2J^;  cos  f„. 

(3) 


Fig  1  Schemalic  0/  infinile  uniforir  plane. wavt  incidtnl  upon  a 

phaae  acreen 
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rociprocoJ  paths,  S ,<  ■  I2:,  ■  0.  and  when  w«  average 
over  many  realisations  of  the  phase  screen  Ekj  (3)  reduces  to 

(/)  -  4/0  4  2/0  +  2/0  +  2/0  <coe  »,j)  4  2/0  <co*  »,j>. 

</)  -  4/0  4  (2/o  4  2;,.),  (4) 

</)  *  inooberent  portion  4  coherent  portion,  (5) 

where  (cos  #12)  •  <ooe  (ij)  -  0  because,  for  directions  other 
than  that  of  the  incident  wave,  all  phase  differences  are 
random. 

Recently  Jakeman  theoretically  modeled  the  scattering 
that  occurs  to  an  optical-beam  wave  on  a  double-paas  path 
through  the  same  phase  screen.  In  bis  analysis  he  assumed 
the  following; 

(1)  The  induced  spatial-phase  fluctuation  e  by  the 
screen  is  a  spatial-Gauaaian  process  with  a  smoothly  varying 
spatial  spectrum. 

(2)  The  mean-square  phase  fluctuation  4,1^  is  much 
greater  than  unity,  »  1,  so  that  the  scattering  is  strong. 

(3)  The  phase  screen  is  illuminated  by  a  monochromatic 
Gaussian- beam  wave  with  radius  which  is  taken  to  be 
much  larger  than  the  largest  diameter  of  index  inhomogene- 
ily.  Hence  the  scattering  gives  rise  to  a  complex  Gaussian 
spatial  process  in  the  far  field. 

(4)  The  plane-reflecting  mirror  is  li'ccted  within  the 
Fresnel  region  where  the  illuminated  area  is  much  greater 
than  the  first  Fresnel  zone. 

(5)  The  plane  mirror  is  located  at  a  distance  z  beyond  the 
phase  screen  within  the  Fresnel  region  of  the  beam  .such  that 
2  «  4H‘i,V4,  where  k  »  2»/A. 


-  /,"'exp^-^  j  |1  4  exp(-2>r'/-Q-)''-),  (9) 

where  P  reduces  to  (drrvi)’.  The  leading  exponential  factor 
daacribes  the  intenaity  angular  diatribution  of  the  acatiered 
waves  oaring  to  two  independent  phase  screens  with  identi¬ 
cal  spatial  apectra.  The  exponential  term  arithin  the  brack - 
eta  deacribes  the  intensity-enhancement  distribution. 

4.  SIMULATION  ALGORITHM 

The  wave  field  is  expressed  by 

£(M  “ /(Mexp(-/*M,  (I0> 

where  f(r)  is  the  complex  amplitude,  k  «  2»/X  is  the  free- 
tpace  wave  number,  and  r  ■  (i.y,  z)  is  the  radius  vector.  If 
one  assumes  propagation  in  the  r  direction,  the  parabolic 
wave  equation  becomes 

2>*  ^  /(M  4  V,V(r)  4  2k■n,{t^f^^,t)  -  0.  (Ill 

02 

where  Oj  *  n  -  (n)  is  the  index  of  refraction,  and  v,-  «  (d-/ 
dx’l  4  (d^/dy’l  is  the  transverse  Laplacisn  In  free  space, 
Eq  (111  assumes  the  form 

2j* -r /(.<•  r)  ■ 'v, -/(.«,  r).  (I2i 

az 

where  .<  »  (i,  y)  is  the  two-dimensional  spatial  variable 
The  Solution  of  the  modified  wave  equation  is  found  liv 
performing  a  Fourier  transform  on  both  sides,  thereby  re¬ 
ducing  Eiq.  ( 121  to  an  ordinary  differential  equation 


Using  these  assumptions.  Jakeman  derived  an  expression 
for  the  far-field  intensity  distribution  of  the  backscattered 
wave.  He  showed  that  there  is  an  enhancement  of  the  back- 
acattered-lighl  intensity  in  the  direction  of  the  incident 
beam.  For  a  mirror  placed  at  a  disunce  sufficiently  far  from 
the  screen.  Gaussian  speckle  will  be  developed  at  the  mirror. 
It  was  shown,  for  this  case,  that  the  average  angularly  dis¬ 
tributed  interuily  for  a  Gaussian  beam  is 


ne.) 


1  4 . 


W',’ 


tip(-2k-p^)''^ 


(61 


where  is  the  peak  inleruity  of  the  backscattered  beam,  d 
is  the  distance  from  the  screen  to  the  mirror,  and  rrtf,  repre¬ 
sents  the  rms  value  of  the  sine  of  the  scattering  angle  of  the 
phase  screen.'  Then  the  term  <f-m„-  is  proportional  to  the 
area  of  the  scattering  disk  at  the  screen  on  the  second  pass. 
The  other  parameters  are  given  by 

+ (dm., (7) 

G  “  '^/jdin  ain  P,  .  P  «  V,(sin  P.  4  sin  P,/,  (8) 

where  P,  is  the  observed  sc*'tiering  angle  and  P,  is  the  angle  of 
incidence. 

A  reasonable  approximation  to  Eq.  (6J  for  the  case  of  a 
plane  wave  would  be  M'„  »  dm.,  so  that 


dF(t.  z) 
dz 


2k 


FU.O). 


(13) 


with  spatial  frequency*  “  (*,.*, I, and  the  .solution  given  bv 

£(»,  z)  *  £(*,  01exp(j*'z/'2A),  (U| 

where  expO« -z /2k )  represents  the  F resnel  factor.  Using  Eq . 
(14)  and  lU  inverse,  we  may  transform  from  the  frequency 
domain  to  the  apace  domain  to  as  to  simplify  calculations  of 
phase-screen  scattering. 

A  two-dimensional  random  phase  screen  is  generated  that 
consists  of  Af  X  Af  equiapaced  real-phase  shifu  «(i),  such 
that  the  screen  is  a  member  of  an  ensemble  having  a  speci¬ 
fied  expected  phase  spectrum  t,!*).  The  average  power 
spectrum  of  the  phase  fluctuations  from  an  ensemble  of 
such  screens  converges  to  a  apecified  function. 

The  phase  correlation  function  of  the  screen  is 

RI4)  ■  1 1  (n,(r.*|ln,(z',l2))  dzdz',  (151 

R(t)  ■  k‘AzA„(s,  -  (161 

where  Az  is  the  width  or  thickness  of  the  screen  A,,  repre¬ 
sents  the  correlation  function  and  may  be  expressed  in  terms 
of  the  phase  specirum'- 

A(*)  “  2»  I  j  4’„(*_.  •«  0, » lexpii*' ■  .< Id'*.  (ITi 

Fourier  transforming  Eq.  (16)  and  using  Eq  (ITl.  we  find  the 
spectrum  of  the  phase  Uuciuations  to  be 
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♦,(«)  -  2*k=A,».(..  -  0.  •).  (18) 

The  diecrete  phaae  screen  is  described  by  s  recUnfle  with 
■idss  L,  •  Mtx  and  L,  ■  NAy,  where  Ax  and  A^  are  the 
spatial -aaiBpling  intervals  on  the  screen.  (For  our  study 
they  ate  taken  to  be  equal.)  Pseudorandom  A(«)  are  gener- 
ated  to  satisfy 

-  LA ♦.<*). 

where  #(«)  is  the  transformed  phase  and  4>«(<)  is  the  phase 
power  spectrum.  By  generating  random  phase  ahifu  that 
satisfy  Eq.  (19)  and  by  taking  the  two-dimensional  discrete 
Fourier  transform  {2DFFT),  we  And  that  the  result  is  a 
phase  screen  that  is  s  desired  member  of  the  ensemble.  In 
summary,  the  phase  screens  are  created  by  filtering  white 
Gauaaian  noiae  to  obtain  a  random  field  erith  the  desired 
spectral  properties  jrandom  phase  shift  p  ■  (it  -f 
where  A  and  B  are  independent  random  variable*  that  are 
realisations  of  a  normally  distributed  random  variable  that 
has  a  variance  e-  •  L^L, <',(•)].  A  power-law  spectrum  that 
contains  power  at  all  scales  was  used.  Any  discrete  simula¬ 
tion  can  contain  only  a  finite  range  of  scales.  Thus  there  are 
inner  and  outer  scales  implicit  in  any  discrete  formulation 
owing  to  the  finite  number  of  sample  points  and  the  finite 
sampling  intervals. 

The  next  concern  is  creating  a  mirror  so  that  the  propagat  ■ 
ing  wave  may  be  reflected  back  through  the  same  phase 
screen.  This  is  accomplished  through  a  window  acting  as  a 
filter  that  represents  the  aperture  describing  a  mirror.  The 
mirror  was  chosen  to  be  so.me  uiukiple  of  the  sampling  inter¬ 
val  of  the  screen.  Calculations  showed  that,  as  long  as  the 
mirror  siie  was  greater  than  6.38Aj  X  6.38A.V,  we  would  be  in 
the  Fresnel  region  of  the  mirror. 

We  summarite  with  the  following: 

(1)  Create  random  phase  shifts  «($ I. 

(2)  Perform  the  2DFFT  on  exD[/tf(.*>i].  which  gives  the 
transformed  field.  £(>,  0). 
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(3)  Multiply  by  the  Fresnel  factor,  exp(;«^*/2*).  for 
propagation. 

(4)  Reflect  the  complex  field  from  the  mirror. 

(5)  Propagate  back  to  the  original  phase  screen  by  multi¬ 
plying  by  the  Fresnel  factor. 

(6)  Perform  the  inve-se  2DFFT  to  obtain  the  field  Eu. 

*). 

(7)  Multiply  £(f,  r)  by  expLrXS)]. 

(8)  Perform  two-dimensional  discrete  Fourier  transfor¬ 
mation  to  obtain  far-field  intensity. 

The  above  algorithm  is  graphically  ahowm  in  Fig.  2.  This 
algorithm  leads  to  the  study  of  a  wave  reflected  in  the  back¬ 
ward  direction.  In  order  to  see  the  enhancement,  one  must 
study  the  intensity  statistics  in  the  far  field,  which  implies 
that  the  Fourier  transform  of  the  electric  field  be  taken  at 
the  exit  plane  of  the  second  screen  and  averaged  on  an 
intensity  basis. 

The  simulation  has  been  developed  so  that  all  propagation 
lengths  are  normalized  with  respect  to  the  sampling  interval 
A,.  Let  M  be  the  numlier  of  samples  such  that  the  spatial- 
frequency  resolution  A«,  is  given  by  the  expression 

As  -  2vM/A,.  (20) 

The  power-law  spectrum  for  the  phase  screen  is  a  function 
of  the  strength  of  turbulence  denoted  by  the  parameter  T. 

4',(i)  -  2»T.".  (21) 

where  n  is  the  power-spectral  index.  The  sirenplh  T  mav  lie 
written  as  a  function  of  U,  which  is  taken  to  be  the  square  of 
the  sciniillaiion  index: 

L;--y,(l-''o/2)sin(i^)(^J"  '. 

An  understanding  of  the  boundaries  involved  for  the  screen 
is  obtained  from  the  mutual  coherence  function: 

R(J)  -  eipl-VA(*)l.  <22) 


Fig  2.  Schematic  of  program  simulation  and  cnmputauon  algorithm 
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where  Dit),  the  wave  structure  function  for  power-law  form, 
is  fivtn  by  the  following  expression: 


which  was  derived  from  the  solution  to  the  fourth  moment  of 
Eq.  Solving  for  R(t)  gives  DU)  results  in 


RU)  “  expl 


1 

2(a 

-  2)r(^ 

:Wt) 

(241 


Hence  the  larger  the  U,  the  narrower  RU)  becomes  and  the 
more  finely  sampled  the  screen  must  be.  Solving  for  the 
angular  spectrum  of  plane  waves  by  inverse  Fourier  trans¬ 
forming  RU),  we  are  able  to  obtain  the  parameters  needed 
for  the  simulation.  These  include  the  sampling  iniervai 
taken  to  be 


A>-vXi/16.  (26) 

Hence  Ly  *  I6(Xr)''7’,  and  we  have  Lx  *  2L\.  where  the 
chosen  screen  size  i*  512  x  266  with  an  1 1 /.I  spectral  index 
The  above  parameters  have  been  validated  and  compared 
with  other  screen  sizes  by  Filice.  The  parameters  that  we 
held  ainslani  for  our  experiment  are  the  aperture  width, 
which  is  taken  to  Ije  a  KKJAr  by  lOOiy  square  centered 
within  the  obserx'inp  area;  the  wavelength,  which  is  taken  to 
be  1  *im;  and  the  distance  set  Ui  20  m.  which  is  past  the  focus 
disunce  of  the  screen. 

The  testing  of  the  program  was  conducted  by  using  a 
deterministic  phase  shift  Ci  >  0.  which  resulted  in  a  sinc- 
funciiun  far-field  pattern.  When  a  linear  phase  is  inserted, 
the  result  is.  across  the  screen,  a  sine  function  shifted  by  the 
known  amount. 


5.  RESULTS 

The  theory  of  enhanced  backscatter  relies  on  the  wave's 
passing  through  the  original  screen  that  it  traversed,  result¬ 
ing  in  a  coherent  addition  of  fields.  On  the  other  hand  when 
the  phaae  acreeru  are  independent,  we  have  an  incoherent 
addition  and  therefore  no  enhancement.  In  Fig.  3  the  algo¬ 
rithm  was  altered  so  as  to  produce  two  independent  phase 
acreeru.  The  parameters  for  this  case  are  80  realizstioru 
with  Lf  w  10.  As  predicted,  there  is  only  normal  scattering 
vrith  ito  enhancement,  thereby  producing  another  roearu  of 
validating  the  program.  The  next  run  coruists  of  two  identi¬ 
cal  acreeru  to  model  the  backscatter.  This  cose,  shown  in 
Fig.  4,  duplays  the  enhancement,  which  it  a  factor  of  2  above 
the  normal  for  V  "  10, 1(X3  realizations.  The  same  case  was 
run  for  V  •  5.  shown  in  Fig.  5.  At  expected,  by  lowering  the 
strength  of  turbulence,  we  widened  the  enhancement  width 
The  last  study  conducted  for  this  paper  involves  the  angle  of 
incidence  of  the  incoming  wave.  VVe  altered  the  incident 
angle  by  0.012  deg.  which  resulted  in  the  shifting  of  the 
enhancement  by  the  tame  amount  shown  in  Fig.  6.  This  run 
was  for  110  realizations  with  (7*10. 

The  cases  that  involve  independent  screens  and  strength 
of  turbulence  U  ■  5  were  run  on  a  Gould  Concepts  32'f7 
computer  with  an  8600  image-processing  board  The  ex¬ 


periments  took  76  and  57  h,  to  run  for  independent  screens 
and  U  w  5,  reapectively.  The  study  of  changing  the  angle  of 
incidence  was  conducted  on  a  DEC  VAX  computer  and  took 
77  h.  The  lost  case,  for  U  *  10  at  100  realizatioru,  waa 
performed  on  the  San  Diego  Cray  XMP/48  computer  and 
took  '>'10  min. 

6.  DISCUSSION 

At  yet  there  does  not  exist  a  tractable  analytical  formulation 
of  the  backscattered  wave  problem  for  a  phase  screen  with  a 
p«iwer-law  spatial  spectrum  of  index  inhumugenieties.  Al¬ 
though  Jakeman’s  analysis  dealt  with  the  scattering  phe¬ 
nomenon  due  to  a  Gaussian  spatial-correlation  function,  it 
can  be  considered  to  be  an  approximation  to  a  power-law 
spectrum  with  a  finite  outer  scale,  i.e.  a  lowest  spatial  fre¬ 
quency.  The  simulation  performed  here  included  an  inner 
and  an  outer  scale  only  through  the  maximum  grid  size  and 
minimum  sparing  between  pixels.  The  Gaussian  power 
apectrum  is  defined  by  a  single  inhnmogeneity  scale,  so  if  the 
illumination  of  the  random  phase  screen  greatly  exceeds  this 
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Fig  3.  Angular  disuibulic  i  of  the  far-field  iniensiiy  of  the  bacli. 
scatter  for  two  independenl  phase  screens,  with  (.'  •  10.  X  -  1  sm.  d 
w  20  m.  and  mirror  size  ■  lOOAz  X  lOOAy  for  80  resiizalions 
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Ftg  4.  Anrulsr  ditlribuuon  of  the  (sr  field  inveniii>  of  the  bscli 
scatter  for  s  normally  incident  wave  with  I'  •  10.  A  •  1  um.  c  ■  20  m. 
and  mirror  size  ■  ](XIAx  X  lOOAy  for  1(X)  realizauonz 
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Fif .  5.  AafuUr  diatribution  of  the  far-rield  inUMity  with  the  tame 
parameuw  (ivan  in  Fi(.  4  for  U  “6. 


4  exceot  (or  110  rcaliutiotu 


scale  the  scattered  optical  field  can  truly  approach  Gaussian 
speckle.  The  scattering,  which  is  due  to  a  power-law  spatial 
spectrum,  approaches  a  Gaussian  speckle  as  a  limit  distribu- 
uon  only  at  an  infinite  distance  from  the  screen  because  of 
the  presence  of  infinitely  large  scales  of  inhomogeneities. 
The  Gaussian  power-law  spectrum  also  may  manifest  strong 
caustics,  whereas  a  true  power-law  spectrum  does  not  do  so. 

The  eridth  of  the  enhancement  peak  for  the  weaker  screen 
is  approsimately  twice  that  of  the  strong  screen.  This  pro¬ 
portionality  of  the  width  to  the  scattering  strength  is  coiuit- 
tent  arith  the  prediction  for  the  Gaussian  spectrum  screen  as 
given  by  Eq.  (9).  The  scattering  disk  is  the  area  of  the  phase 
screen  that  encompasses  all  the  scattering  that  contributes 
to  a  particular  point  at  the  mirror.  The  radius  of  this  disk  is 
given  by  6^,  where  0,  is  the  sridth  of  the  angular  spectrum 
and  r  is  the  distance  to  the  mirror.  The  angle  0,  •  l/kso. 
where  k  «  2vA  and  s<.  is  the  width  of  the  mutual  coherence 
function  defined  by  equating  the  plsne-wave  structure  func¬ 
tion  to  2.  i.e..  £)p(sti)  “  2.  For  our  case  then  we  have 


The  angular  sridth  of  the  inteiuity-enhancement  pattern  is 
deUrminad  by  the  diffraction  width  of  the  scattering  disk. 
Only  rays  separated  by,  at  most,  the  diameur  of  the  scatur- 
ing  disk  can  have  reciprocal  paths  and  add  coherently. 
Hence  the  scattering  disk  determines  the  equivalent  aper¬ 
ture  for  the  enhancement  and  thereby  its  angular  width. 
The  far-fwld  enhancement  width  is  0,  >  2.44  X/D,  where  D  is 
the  scattering  disk  diameter,  D  •  20^.  Theoretically,  there 
is  s  factor  of  1.5  from  10  toU^Siti  the  width  of  the 
enhancement.  Measuring  the  width  at  50%  below  the  peak 
so  as  not  to  have  contamination  from  the  normal  scattering, 
we  find  that  the  factor  between  the  two  strengths  is  -^2.0 
This  difference  may  be  contributed  to  the  additional  diffrac¬ 
tion  from  the  finite  mirror  aize  of  the  simulation.  Cslrula- 
tions  show  that  the  diffraction  from  the  mirror  is  -^30%  of 
the  theoretical  enhancement  width.  Hence,  by  taking  this 
factor  into  account,  we  find  that  the  width  agrees  with  the 
results  from  the  numerical  experiment. 

For  both  the  weak-  and  strong-scattering  screens  the 
height  of  the  enhancement  peak  is  approximately  twice  that 
of  normal  scatter.  The  enhancement  predicted  by  Eiq.  (9i  in 
the  backward  direction,  Q  •  0,  is  twice  that  of  the  incoher¬ 
ently  scattered  angular  pattern.  There  seems  to  be  good 
agreement  between  the  theory  for  the  Gaussian  spectrum 
screen  and  the  power-law  screen  of  the  simulation  Hence 
the  essence  of  the  phenomenon  does  not  seem  to  be  depen¬ 
dent  on  the  details  of  the  phase-screen  spectrum. 
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5.  IR  FOCAL  PLANE  ARRAY  TESTING 


Glenn  Boreman,  Ph.D 

Infrared  Systems  Laboratory /CREOL 

5.1  Background 

CREOL’s  Focal  Plane  Array  testing  program  consists  of  the  measurement  of  various  specific 
detectors,  the  development  of  various  instrumentation  devices  for  the  measurement  of  detector 
modulation  transfer  functions  (VfTFs),  the  development  of  an  IR  scene  projector  system,  as  well 
as  the  development  of  data  processing  and  computing  facilities. 

T^e  types  of  detectors  that  have  been  tested  include  an  ISTEF  provided  Mitsubishi  Platinum 
Siiicide  CCD  camera,  a  CREOL  Mitsubishi  PtSi  CCD,  a  Signal  Processing  In  The  Element 
(SPRITE)  Detector,  a  silicon  CCD  staring  array,  and  hexagonal  array  structures. 

Unique  instrumentation  has  been  developed  for  the  measurement  of  the  MTF  of  three 
detector  architectures:  1)  self-scanned  arrays,  2)  SPRITE  detectors,  and  3)  staring  arrays. 

An  IR  scene  projector  instrument  has  been  developed  to  provide  complex  realistic  images 
to  the  system  array  under  test.  The  goal  of  this  system  is  to  provide  a  realistic  time  history, 
dynamic  range,  and  resolution  in  response  to  an  externally  generated  video  source. 

The  computing  facilities  have  been  upgraded  from  an  IBM  PC  AT  to  a  SUN  330 
workstation  with  future  plans  to  implement  a  fiber  link  to  the  University’s  SUN  network. 

5.2  Characterization  of  Mitsubishi  Platinum  Siiicide  CCD  Camera 

An  important  component  of  the  infrared  imaging  system  at  ISTEF  is  the  512x512  PtSi 
charge-coupled  device  (CCD)  array  from  Mitsubishi.  This  camera  was  borrowed  from  ISTEF  on 
two  occasions.  Extensive  test  data  were  collected  during  that  limited  period.  TTie  data  from  those 
measurements  yielded  calibration  information  for  the  camera.  Specific  tests  performed  measured 
uniformity,  linearity,  responsivity,  and  mean-variance. 

The  experimental  setup  for  each  characterization  used  a  2.5-mm-diameter  1000  K  blackbody 
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source.  The  distance  from  the  blackbody  source  to  the  CCD  array  was  varied  to  produce  a 
complete  range  of  irradiances  on  the  array.  The  small  source  size  produced  uniform  illumination 
over  the  detector  array  (flat-field  inputs).  No  lens  was  used  in  this  setup;  thus  the  flux  from  the 
blackbody  followed  a  well-defined  inverse-square  dependence.  The  passband  of  the  system  was  3 
to  5  ^nt,  defined  by  the  Dewar  window. 

5.2.1  Uniformity 

Figure  5-1  is  an  example  of  the  raw  array  output,  with  a  flat-field  input.  The  spatial  noise 
pattern  is  typical  of  uncorrected  output  from  PtSi  CCDs'  and  represents  a  3%  nonuniformity 
(standard  deviation  divided  by  maximum  dynamic  range).  When  the  internal  camera  calibration 
is  enabled  (correction  for  gain  and  offset),  the  output  produced  for  a  flat-field  input  is  uniform  to 
well  below 

5.2.2  Linearity 

Figure  5-2  illustrates  response  in  pixel  count  plotted  against  focal-plane  irradiance  in  W/cm‘. 
The  camera  response  is  linear  at  high  irradiance  levels,  until  the  saturation  limit  of  the  camera  is 
reached.  Minor  departures  from  linearity  occur  at  low  irradiance  levels,  and  can  be  corrected  by 
using  look-up  tables. 

5.2.3  Responsivitv 

The  slope  of  the  curve  in  the  linear  regime  yields  a  responsivity  of  1.5  x  lO  ’  W  cm'^  pixel 
count  ',  which  agrees  with  the  Mitsubishi  specification  for  responsivity. 

5.2.4  Mean  Variance 

For  each  frame  of  data  representing  a  particular  blackbody  distance,  the  variance  of  the 
pixel  count  is  plotted  against  the  mean.  Figure  5-3  shows  the  resulting  mean-variance  curve, ^ 
which  indicates  whether  or  not  the  array  is  operating  under  background-limited  (BLIP) 
conditions.  If  the  primary  contribution  to  the  variance  is  the  photon  noise  on  the  background, 
then  the  plot  of  mean  versus  variance  will  be  linear,  consistent  with  Poisson  statistics.  Our  data 
in  Figure  5-2  indicate  that  the  Mitsubishi  CCD  imager  is  not  BLIP  for  focal-plane  irradiances 
below  approximately  20x1  O'*  W/cm^. 

To  properly  support  this  program,  we  purchased  a  Mitsubishi  PtSi  array  identical  to  that 
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residing  at  ISTBF.  Thus  we  can  continue  the  characterization  of  f>erformance  without 
interruption  of  ISTEF  data-acquisition  activities.  The  acquisition  of  this  equipment  was 
supported  by  CREOL  research  initiation  funds. 

53  Focal  Plane  MTF  Research 

The  modulation  transfer  function  (MTF)  is  perhaps  the  most  useful  figure  of  merit  for 
electro<optical  systems.  A  ^>ecification  of  image  quality  in  the  spatial-frequency  domain 
provides  a  flexibility  and  an  ease  of  analysis  unmatched  by  a  spatial-domain  description.  MTF 
allows  the  convenient  cascading  of  various  subsystems  by  a  simple  multiplication  of  their  transfer 
funaions,  rather  than  by  the  convolution  of  an  equivalent  number  of  impulse  responses. 

Focal  plane  arrays  remain  one  important  subsystem  for  which  the  application  of  MTF 
techniques  is  still  under  development.  While  classical  optical  systems  are  well  described  by 
conventional  MTF  theories,  focal  plane  arrays  have  a  number  of  complicating  factors  related  to 
aliasing,  finite  sampling  windows,  charge  transfer  inefficiency,  crosstalk,  and  system  noise.  In 
addition,  even  the  instrumentation  necessary  for  the  measurement  of  MTF  for  focal  plane  arrays 
is  still  a  subject  of  research. 

The  research  on  MTF  of  focal  plane  arrays  at  CREOL  separates  into  three  areas,  each 
of  which  will  be  discussed  separately.  These  include;  development  of  MTF  instrumentation  for 
self-scanned  arrays;  development  of  M  J  h  instrumentation  for  staring  arrays;  and  analysis  of  the 
hexagonal-array  MTF. 

5.3.]  MTF  Instrumentation  Development  -  Self-Scanned  Arrays 

SPRITE  detectors  are  self-scanned  detector  structures  that  are  usually  found  in  linear 
focal  plane  arrays.  SPRITE  technology  is  currently  being  used  in  a  number  of  infrared  imaging 
systems  in  the  8-10  waveband  as  an  interim  solution  until  staring  focal-plane  technology 
matures  for  HgCdTe.  SPRITE  technology  represents  a  viable  alternative  for  ISTEFs  long-wave 
infrared  imaging  needs. 

A  SPRITE  detector,  seen  in  Fig.  5-4,  consists  of  a  strip  of  photoconductive  material 
operated  under  a  constant  current  bias.  Photogenerated  carriers  w’ill  drift  along  the  structure  at 
a  rate  determined  by  the  applied  voltage  and  the  carrier  mobility.  If  the  IR  image  is  scanned 
across  the  deteaor  structure  at  the  same  rate  as  the  carrier  drift  velocity,  the  signal-to-noise 
ratio  is  increased  because  of  signal  averaging. 
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Figure  5-4.  SPRl  i  t  detector  configuration 
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During  the  past  year,  we  developed  instrumentation  (Fig.  5*5)  to  measure  the  MTF  of 
SPRITES  by  means  of  the  sine-wave  response.  An  infrared  He-Ne  laser  (wavdenglh  3.39  is 
used  in  omijunction  with  a  set  of  double-slit  apertures  to  project  interference  fringes  of  variable 
^tiat  frequency  onto  the  SPRITE.  The  fringes  are  scanned  at  an  appropriate  velocity  with  a 
polygon  scanner.  The  MTF  is  measured  from  the  modulation  depth  of  the  resulting  signals  as  a 
funaion  of  frequency. 

Our  results  on  the  MTF  of  SPRITE  deteaors  (Fig.  5-6)  have  confirmed  an  expected 
theoretical  dependence’  of  MTF  on  the  length  of  the  detector  (longer  detector  yields  poorer 
MTF  because  of  charge  diffusion).  We  have  submitted  an  article*  on  this  measurement  method, 
which  is  included  as  Appendix  A. 

This  type  of  target  generator  will  also  be  useful  to  the  ISTEF  mission  in  unconventional 
imaging.  The  developmental  system  to  evaluate  the  performance  of  the  unconventional  imaging 
sj'stems  will  require  a  means  to  generate  interference  fringes  of  variable  spatial  frequenc)'  and 
contrast  level. 

5.3.2  MTF  Instrumentation  Development  -  Staring  Arrays 

Staring  focal  plane  atTa>'s  are  a  key  technology'  element  for  the  infrared  imaging  s>'stems 
of  the  future.  The  advantages  in  compactness,  size  and  weight,  mechanical  simplicity,  and 
sensitivity  indicate  that  this  class  of  image  receiver  will  be  the  new  standard,  surpassing  the 
scanned  FLIR  systems  so  prevalent  today.  Indeed,  for  many  ISTEF  missions  b  target  location, 
imaging  spectroscopy,  and  unconventional  imaging  staring  arrays  are  the  central  component  to 
determining  the  ultimate  performance  of  these  systems. 

We  have  developed  a  new  technique  to  measure  the  MTF  of  staring  focal  plane  arrays, 
both  in  the  visible  and  infrared  portions  of  the  spectrum.  The  basic  technique’  consists  if 
projecting  a  laser  speckle  of  known  spatial  frequency  content  onto  the  focal  plane  array.  The 
spatial  power  sp>ectral  densities  (PSDs)  of  the  speckle  patterns,  before  and  after  deteaion,  are 
related  by  the  MTF  of  the  CCD  array: 

where  {  is  the  spatial  frequency'.  We  have  developed*  an  extension  of  the  basic  method  for  speckle 
generation  that  employs  an  integrating  sphere  as  the  phase-randomizing  medium,  rather  than  a 
ground-glass  di^ser.  Figure  5-7  is  a  schematic  of  the  measurement  setup.  The  integrating-spherc 
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Figure  5-5.  MTF  instrumentatioD  developed  for  SPRTTEs 


Figure  5-6.  Results  of  MTF  measurement  for  SPRITES. 
Note  the  decrease  in  MTF  for  longer  SPRITE  element  length 
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approach  is  particularly  useful  at  infrared  wavelengths. 

Our  current  research  in  this  area  involves  the  development  of  techniques  that  allow  the 
MTF  to  be  measured  past  the  ^atial  Nyquist  frequency  of  the  focal  plane  array.  This  requires  a 
capability  to  generate  a  ^leckle  pattern  with  a  flat  bandpass  characteristic  to  avoid  artifaas  in  the 
data  that  would  result  from  aliasing.  Briefly,  the  PSD  of  the  speckle  is  proportional  to  the 
autocorrelation  of  the  aperture  that  generates  the  speckle.’  We  use  a  nonseparable  aperture  (Fig 
5-8)  that  has  a  substantial  flat  region  of  autocorrelation  (Fig.  5-9)  to  produce  the  flat  bandpass.  In 
Fig.  5-10  we  show  the  full  two-dimensional  PSD^,  which  has  been  modified  by  the  CCD  MTF.  We 
are  currently  developing  the  technique  using  a  silicon  CCD  at  a  wavelength  of  0.8  fim,  and  have 
purchased  components  to  extend  these  tests  to  the  platinum  silicide  camera  at  3.39  fim.  An  article 
is  being  prepared  on  this  research  for  submission  to  Optical  Engineering. 

In  addition  to  the  importance  of  laser  speckle  to  the  MTF  research,  the  availability  of  a 
caUbrated  speckle  generator  will  prove  useful  to  the  ISTEF  mission  in  unconventional  imaging.  Our 
techniques  facilitate  the  production  of  speckle  patterns  over  a  wide  range  of  wavelengths  that  have 
well-known  statistics  and  well-defined  spatial-frequency  contents.  One  component  of  the 
developmental  system  to  evaluate  the  performance  of  ISTEF’s  unconventional  imaging  systems  will 
be  a  speckle  generator  using  the  integrating  sphere  configuration. 

5.3.3  Analysis  -  Hexagonal  Staring  Arrays 

We  analyzed  and  compared  MTFs  obtained  from  a  hexagonal  array  of  hexagonal  detectors 
and  a  hexagonal  array  of  rectangular  detectors  (Fig.  5-11).  Hexagonal  detector  arrays  are  valuable 
in  supcr-ccniroiding  applications,*’  facilitating  more  efficient  and  robust  super-centroiding 
algorithms  because  of  the  symmetric  nature  of  the  lattice. 

We  considered  the  spatial  averaging  that  occurs  over  the  surface  of  the  photosite,  and  cast 
it  as  a  convolution  of  the  spatial  photosite  responsivity  and  the  incident  irradiance  distribution. 
Then,  operating  in  the  Fourier  domain,  the  MTF  attributable  to  this  spatial  averaging  is  seen. 

Our  results  are  shown  in  Fig.  5-12.  The  MTF  in  the  vertical  direaion  is  enhanced  (•  10% 
MTF  increase  at  Nyquist),  at  the  expense  of  the  horizontal  MTF  (  •  5%  MTF  decrease  at  Nyquist). 

The  MTF  of  these  structures  has  not  been  previously  presented  in  the  literature,  and  we 
have  submitted  an  article”  xo  Applied  Optics  as  a  result  of  this  research  (see  Appendix  B). 
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Figure  5-7.  Schematic  of  MTF  measurement  setup. 
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Figure  5-8.  Nonseparable  aperture  for  the  sf>eckle  generator 
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Inp'jt  Power  Spectrum  of  Laser  Speck'‘=? 
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Figure  5-9.  Input  power  spearal  density  for  the  laser  speckle  which  results  from  the  use  of 
the  aperture  in  Fig.  5-8.  Note  the  Oat  bandpass  region.  L  is  the  largest  separation  in  the 
aperture,  1  is  the  width  of  the  slits,  Z  is  the  distance  from  the  aperture  to  the  focal  plane,  and 
X  is  the  wavelength. 


Figure  5-10.  Two  dimensional  power  spectral  density  of  the  laser  speckle,  after  processing  by 
the  CCD  array.  The  portions  of  the  spectrum  removed  from  the  central  cross  structure  are 
the  flat  bandpass  portions  of  the  speckle  seen  in  ID  in  Fig.  5-9.  The  shading  seen  in  these 
sections  is  caused  by  the  MTF  of  the  array. 
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Figure  5-11.  Geometry  of  the  hexagonal  arrays  investigated. 
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Figure  5-12.  M  IF  results  for  the  hexagonal  detectors  and  rectangular  detectors  seen  in  Fig. 
5-11.  The  hexagonal  detectors  have  a  larger  MTF  in  the  venical  direction  than  do  the 
rectangular  detectors.  The  rectangular  detectors  have  a  smaller  MTF  advantage  in  the 
horizontal  direction. 
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5.4  Target  Prqjeclor  Development 

A  central  issue  facing  the  users  of  infrared  s)’siems  is  end-io-end  performance 
including  the  optics,  focal  plane  array,  electronics,  and  data -processing  algorithms. 
The  desired  test  images  are  becoming  more  complex.  The  usual  bar  chans,  impulse  response  tests 
and  broad-area  thermal  sensitivity  measurements  must  be  coupled  with  the  pre.seniation  of  realistic 
infrared  scenes  to  the  ^tem  under  lest.  Infrared  scene  projection  is  an  effective  tool  for 
system-performance  diaracterization,  without  the  need  for  expensive  and  time-consuming  field 
trials. 

The  central  technical  problem  is  to  create  infrared  scenes  of  realistic  lime  history,  dynamic 
range,  and  resolution  in  re^nse  to  an  externally  generated  video  source.  We  are  vigorously 
pursuing  research  that  shows  promise  of  filling  this  need. 

The  projector  system  seen  in  Fig.  5-13  (on  long-term  loan  to  CREOL  from  the  Naval 
Training  Systems  Center)  is  a  raster-scanned  COj  laser  operating  at  10.6  >im.  TTie  laser  is 
acousio-optically  modulated  in  response  to  video  rate  input,  and  a  two-dimensional  beam  .scan  is 
performed  to  directly  write  the  infrared  scene  into  the  entrance  aperture  of  a  thermal  imager 
^•stem.  The  projeaor  MTF  is  a  function  of  both  the  modulator  drive  format"  and  the  laser  beam 
profile.^ 

Currently  we  are  performing  MTF  measurements  on  this  system  that  investigate  the 
additional  effect  of  acoustic  wave  decay  in  the  modulator  celL  We  have  presented  a  paper  on  this 
topic  at  the  March  1991  SPIE  meeting  on  »erospacc  Sy-stems  in  Orlando. 

In  addition  to  the  utility  of  this  research  for  imager  characterization,  this  target-projection 
technology  will  also  be  used  in  the  test-bed  system  for  evaluating  the  performance  of  unconventional 
imaging  systems  for  ISTEF. 

SJS  Computer  System  Development 

The  charaaerization  of  focal  plane  arrays  requires  digital  processing  of  high-resolution  data 
sets.  Until  last  year,  we  had  performed  this  processing  on  IBM  AT-class  machines.  The  volume 
of  data  dictated  that  we  progress  to  a  more  capable  data-processing  computer.  With  joint  funding 
from  ONR  and  the  Florida  High  Technology  Council,  we  obtained  a  SUN  330  workstation.  TTie 
typical  processing  time  for  the  512x512  data  sets  has  decreased  by  more  than  an  order  of  magnitude. 
We  are  currently  pursuing  additional  storage  media,  and  access  to  the  UCF  campus  network  by 
means  of  a  Obcr-optic  link  from  the  Research  Park.  Also,  we  have  ordered  a  comprehensive 
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Figure  5-13.  COj  Laser-Based  Scene  Projector 
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daU-prooessing/math  software  package  using  CREOL  research  initiation  funds.  This  system 
promises  to  fill  the  needs  of  the  focaJ-plane-array  research  program  in  terms  of  data  acquisition, 
prooessiiig.  di^ilay,  storage,  and  experimental  control  for  the  foreseeable  future. 

The  SUN  format  has  been  used  at  the  main  ISTEF  sites  and  other  ISTEF  research 
groups  at  UCF.  Adoption  of  this  format  will  increase  the  sharing  of  data  and  software  among  the 
participating  groups. 
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